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ABSTRACT 

Aims. High-resolution Very-Long-Baseline Interferometry observations of relativistic jets are essential to constrain fundamental pa¬ 
rameters of jet formation models. At a distance of 249 Mpc, Cygnus A is a unique target for such studies, being the only Eanaroff-Riley 
Class II radio galaxy for which a detailed sub-parsec scale imaging of the base of both jet and counter-jet can be obtained. Observing 
at millimeter wavelengths unveils those regions which appear self-absorbed at longer wavelengths and enables an extremely sharp 
view towards the nucleus to be obtained. 

Methods. We performed 7 mm Global VLBI observations, achieving ultra-high resolution imaging on scales down to 90 pas. This 
resolution corresponds to a linear scale of only ~400 Schwarzschild radii (for Mbh = 2.5 x 10 ^Mq). We studied the kinematic 
properties of the main emission features of the two-sided flow and probed its transverse structure through a pixel-based analysis. 
Results. We suggest that a fast and a slow layer, with different acceleration gradients, exist in the flow. The extension of the accel¬ 
eration region is large (~10'*Rs)> indicating that the jet is magnetically-driven. The limb brightening of both jet and counter-jet and 
their large opening angles (0j~lO°) strongly favor a spine-sheath structure. In the acceleration zone, the flow has a parabolic shape 
(r oc ^0 55±o 07^ q’jjg acceleration gradients and the collimation profile are consistent with the expectations for a jet in “equilibrium” 
(Lyubarsky 2009), achieved in the presence of a mild gradient of the external pressure {p oc < 2). 

Key words, active-jets - high angular resolution 


1. Introduction 

The physics of extra-galactic relativistic jets is one of the 
most exciting and challenging topics in astrophysics. Consider¬ 
able progress has been made in its understanding over the last 
decades. By interpreting jets as fluid phenomena, their proper¬ 
ties on parsec and kilo-parsec scales are quite well explained by 
applying the laws of relativistic hydrodynamics. However, a fun¬ 
damental problem is still unsolved; how can jets be accelerated 
up to Lorentz factors of tens, and how can they be so sharply 
collimated? The answer to these questions must be intimately 
related to the mechanism of jet launching. The necessity of hav¬ 
ing light jets, able to reach relativistic speeds, gradually led to 
the exclusion of a purely hydrodynamic launching, i.e. exploit¬ 
ing the gas pressure gradient in the accretion disk, as this implies 
too high a mass loading. The most successful models are those 
involving the extraction of rotational energy from the accretion 
disk (Blandford & Payne 1982) or of the black hole itself (Bland- 
ford & Znajek 1977) through the magnetic field. Some numer¬ 
ical simulations (McKinney & Gammie 2004; De Villiers et al. 
2005) have suggested that launching from the accretion disk can 
only produce jets with moderate Lorentz factors (T < 3), while a 
process involving the black hole rotation can more easily achieve 
r~10 (McKinney 2006). However it is not ruled out that the in¬ 
ner disk can launch relativistic jets as well (Vlahakis & Konigl 
2004; Komissarov et al. 2007). The possibility of reaching a cer¬ 
tain terminal Lorentz factor depends strictly on the details of the 
process which converts the initially magnetically dominated jet 
into a kinetically dominated one. In theoretical models describ¬ 


ing the conversion of Poynting flux, the mechanisms of acceler¬ 
ation and collimation of the flow show a strong interplay. The 
magnetic field structure is, in this context, a fundamental pa¬ 
rameter, but certainly not the only one. In jets missing a large 
scale (poloidal) magnetic field, i.e. featuring a purely toroidal 
field (Lynden-Bell 1996), an efficient magnetic self-collimation 
seems unlikely (Spruit et al. 1997), and confinement from an 
external medium appears to be the only possibility. The pres¬ 
sure from a poloidal magnetic field, on the other hand, can ef¬ 
fectively confine the flow. However, an additional contribution 
from the external medium is anyhow necessary for preventing 
the terminal speed and collimation of the flow to be reached on 
unreasonably large scales (Tomimatsu 1994; Beskin et al. 1998). 
In accreting systems like AGN, this medium is likely to be the 
magnetized slow wind ejected by the edges of the accretion disk. 

Exact analytic solutions for these problems are difficult to 
obtain if one takes adequately into account the high degree of 
complexity of the system. While numerical simulations can con¬ 
siderably improve our knowledge, direct observations of the fun¬ 
damental parameters of the flow on small scales are equally de¬ 
cisive. With the development of VLBI techniques at millimeter 
wavelengths and of space-VLBI, the spatial resolution achiev¬ 
able in observations of radio loud AGN has dramatically im¬ 
proved over the last few years, allowing us to zoom into the 
innermost regions of extra-galactic jets. Besides the possibility 
of obtaining explicit measurements of speeds, scale of the ac¬ 
celeration and shape of the flow, high resolution VLBI observa¬ 
tions also allow the internal structure of the flow to be exam- 
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Table 1. Log of observations at 7 mm (43 GHz) and characteristics of the clean maps. Col. 1: Day of observation. Col. 2: Array. VLBA - Very 
Long Baseline Array; GBT - Green Bank Telescope; On - Onsala; Nt - Noto; Eb - Elfelsberg; Ys - Yebes. Col. 3: Correlator. Col. 4: Beam FWHM 
and position angle. Col. 5: Peak flux density. Col. 6: Noise. All values are for untapered data with uniform weighting. 


Date 

Antennas 

Correlator 

Beam EWHM [mas, deg] 

Speak[mJy/beam] 

rms [mJy/beam] 

23/10/2007 

VLBA*, GBT, On, Eb, Ys 

JIVE MKIV 

0.23 X 0.11,-21.9° 

197 

0.13 

16/10/2008 

VLBA*, GBT, On, Nt, Eb, Ys 

JIVE MKIV 

0.22x0.10,-11.0° 

264 

0.11 

19/03/2009 

VLBA, GBT, On, Nt, Eb, Ys 

Bonn MKIV 

0.23 X 0.09,-12.5° 

269 

0.14 

11/11/2009 

VLBA, GBT, On, Nt, Eb 

Bonn MKIV 

0.19x0.10,-20.1° 

226 

0.32 


*Saint Croix data not available. 


ined. Recent parsec and sub-parsec scale imaging has revealed 
that relativistic jets are far from the simple assumption of ho¬ 
mogeneous conical outflows, and show instead a complicated 
transverse structure. Among other characteristics, limb bright¬ 
ening was observed both in some nearby FRI radiogalaxies, like 
M 87 (Kovalev et al. 2007) and 3C 84 (Nagai et al. 2014), and 
in blazars, as in the case of Mrk 501 (Giroletti et al. 2004) 
and 3C273 (Lobanov & Zensus 2001). Limb brightening has 
often been interpreted as the result of a spine-sheath structure 
of the jet, i.e. made of a fast and light ultra-relativistic spine 
embedded in a cylindrical flow of slower and denser material 
(Sol et al. 1989). Such a scenario provides a good explana¬ 
tion for some inconsistencies between observations and the stan¬ 
dard unified model for AGN (Barthel 1989), as for example the 
lower Lorentz factors and lower intrinsic jet powers derived in 
FRI galaxies (sheath-dominated) with respect to BL Lac ob¬ 
jects (spine-dominated) (Hovatta et al. 2009; Xu et al. 2000). 
Also, a spine-sheath interaction seems to be an essential ingre¬ 
dient for Spectral Energy Distribution (SED) modeling of TeV 
BL Lacs, where the high energy emission can hardly be recon¬ 
ciled with the low Lorentz factors deduced from VLBI studies 
(Ghisellini et al. 2005). The origin of such stratification is how¬ 
ever unclear. A two-component jet with the described properties 
could directly arise (Hardee et al. 2007; Xie et al. 2012) from a 
jet launching process in which the external sheath is ejected from 
the accretion disk, while the central spine is fueled from the er- 
gosphere. Alternatively, as shown in Lobanov & Zensus (2001) 
for 3C273, a double-rail structure can also be effectively cre¬ 
ated by helical Kelvin-Helmholtz instabilities. In this case, if the 
emission is dominated by the pressure-enhanced regions of the 
instability, the observed motions will be those of the instability 
patterns, typically propagating at much lower speeds (0.1-0.5 c) 
compared to the flow or the shocks in the jet. 

In the following we present a study of the jet transverse 
structure and kinematics of the Fanaroff-Ryley II (FRII) ra¬ 
dio galaxy Cygnus A, from Global VLBI observations at 43 
GHz. At the source redshift (z=0.056, 1 mas~ 1.084 pc, assum¬ 
ing a ACDM cosmology with Ho= 70.5 h“^ km s Mpc 
Qm = 0.27, Q.\ - 0.73), these provide an angular resolution 
of ~400 Schwarzschild radii for Mbh = 2.5 x IO^Mq (Tad- 
hunter et al. 2003), and allow the immediate surroundings of the 
central engine to be imaged. In addition to the proximity, the 
choice of this target is motivated by further advantages. Firstly, 
the jet of Cygnus A is seen at a large viewing angle, in the range 
50° - 85° (Bartel et al. 1995, and references therein). This means 
that both geometrical and relativistic effects are much reduced, 
implying respectively that the intrinsic morphology and speed of 
the flow can be more easily inferred and that a counter-jet can 
be detected. Moreover, both the main jet and the counter-jet can 
be transversally resolved on sub-parsec scales, which is funda¬ 
mental for studying the collimation. Ultimately, while more such 
studies have been carried out on FRI radiogalaxies (e.g M87, 


3C84, CenA), few observational constraints have been provided 
for strong FR II objects, mainly because their jets are generally 
fainter on parsec scales. Cygnus A is in this respect a unique 
target. The paper is organized as follows. In Sec. 2 we present 
the observational setup and data reduction; in Sec. 3 we describe 
the methods used for the kinematic analysis and the study of the 
transverse intensity profiles; in Sec. 4 & 5 we present and discuss 
the results, which are summarized in Sec. 6. 


2. Observations and data reduction 

The data set comprises 4 epochs at 7 mm (43 GHz) from Global 
VLBI observations made between 2007 and 2009 (Tab. 1). The 
long total on-source time of ~ 8 hours, spanned over a ~16 hour 
track, and the large number of antennas (varying between 13 and 
15) produced a good uv-coverage. During observations, antennas 
were frequently switched between the target and the nearby cali¬ 
brator 2013-1-370, in order to reduce pointing errors at the biggest 
dishes. For the first three epochs, data were recorded in dual po¬ 
larization mode using 4 subbands (IFs) with a total bandwidth 
of 64 MHz per polarization (512 Mb/s recording rate). The last 
epoch, instead, was observed in single polarization, with 16 sub¬ 
bands and a total bandwidth of 128 MHz. Data were calibrated 
in AIPS (Astronomical Image Processing System) following the 
standard procedures, including fringe fitting, opacity correction 
for atmospheric attenuation and bandpass calibration. The imag¬ 
ing and self-calibration of amplitude and phase were performed 
in DIFMAP (Difference Mapping). 


3. Data analysis 

3.1. Model fitting and aiignment of the images 

Since the absolute position information is lost during phase self¬ 
calibration of VLBI data, one of the main difficulties in multi¬ 
epoch studies is image registration. This can be particularly chal¬ 
lenging in the case of radio galaxies because, unlike in blazars, 
a prominent VLBI core to be taken as the stationary reference 
point is usually missing. In order to obtain a simplified model 
of the source, we have fitted circular Gaussian components to 
the complex visibilities using the MODELFIT subroutine in 
DIFMAP. The components were then cross-identified in the 4 
epochs, taking into account their size, flux density and the mor¬ 
phology of the region (Fig. 1 & Tab. 2). The images were Anally 
shifted and aligned on the position of component N, assumed 
to be stationary. This feature, although it may not coincide with 
the true jet base, is the most compact and describes a region with 
negligible morphological changes. The same choice was made in 
Krichbaum et al. (1998) (there named JO), based also on a spec¬ 
tral analysis. Note that, by applying the aforementioned shift, 
the positions of a sharp gap of emission at ~ 0.2 mas to the east 
of component N also become aligned. The nature of this feature. 
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already noticed in Bach et al. (2008) is discussed in Sec. 5.6. Ac¬ 
cording to the cross-identification described, a new component 
(J1) was ejected between October 2007 and March 2009, accom¬ 
panied by a brightening by about 50% of the core and inner-jet 
region. The kinematics resulting from this alignment will be pre¬ 
sented in Sect. 4.1. 

Concerning the errors associated with the each parameter of 
the gaussian components, reported in Tab. 2, a rigorous method 
taking into account the local signal to noise in the image around 
each feature was applied (see Lobanov 2005; Schinzel et al. 
2012, for the details of the method). Since the positional errors 
obtained are quite small compared to the beam size, conservative 
values equal to lOcr for the radial separation (Col. 3) and 5cr for 
the transverse size (Col. 5) were assumed. 




Fig. 1. 43 GHz modelfit maps of CygA. Images were tapered and con¬ 
volved with a common beam with FWHM of 0.2 x 0.1 mas, pa=-15°. 
Contours represent isophotes at 1.8, 3.6, 7.2, 14.4, 28.8, 57.6, 115.2 
mjy/beam. The reference point for the alignment is component N, as¬ 
sumed to be stationary. 


jets, first fitting a single Gaussian and then fitting two Gaussian 
profiles. 

Table 2. Properties of the MODELFIT gaussian components from the 
first to the last epoch (1: 2007.81, 2: 2008.79, 3: 2009.21, 4: 2009.86). 
Col. 1: Component label. Col. 2: Epoch. Col. 3: Radial separation from 
the core N. Col. 4: Integrated flux density. Col. 5: Transverse size 
(FWHM). Col. 6: Apparent speed in units of c. 


ID 

Ep. 

z [mas] 

5. [Jy] 

r [mas] 


Cl 

1 

-0.504±0.043 

68±14 

0.240±0.043 



2 

-0.545±0.063 

61±16 

0.220±0.063 



3 

-0.529±0.041 

78±15 

0.243±0.041 



4 

-0.508±0.046 

79±17 

0.222±0.046 

0.00±0.06 

N 

1 

0 

254±20 

0.087±0.002 



2 

0 

223±24 

0.076±0.003 



3 

0 

313±31 

0.101±0.003 



4 

0 

212±14 

0.090±0.002 


JO 

1 

0.131±0.002 

212±15 

0.124±0.002 



2 

0.130±0.003 

473±39 

0.125±0.003 



3 

0.098±0.004 

201±25 

0.083±0.004 



4 

0.118±0.003 

341±27 

0.117±0.003 

-0.03±0.03 

J1 

1 

2 






3 

0.226±0.005 

446±40 

0.151±0.005 



4 

0.357±0.026 

359±55 

0.228±0.026 

0.75±0.04 

J2 

1 

0.278±0.004 

53±7 

0.060±0.004 



2 

0.342±0.028 

44±11 

0.066±0.028 



3 

0.429±0.039 

31±9 

0.060±0.039 



4 

0.621±0.013 

34±6 

0.062±0.013 

0.96±0.05 

J3 

1 

0.385±0.005 

212±17 

0.199±0.005 



2 

0.467±0.040 

143±27 

0.257±0.040 



3 

0.629±0.026 

140±21 

0.252±0.026 



4 

0.789±0.014 

84±13 

0.138±0.014 

L15±0.04 

J4 

1 

0.731±0.051 

108±21 

0.310±0.051 



2 

0.805±0.121 

92±26 

0.332±0.121 



3 

0.896±0.048 

90±17 

0.291±0.048 



4 

0.907±0.078 

94±22 

0.395±0.078 

0.37±0.07 

J5 

1 

1.129±0.113 

50±12 

0.411±0.113 



2 

1.230±0.245 

59±21 

0.432±0.245 



3 

1.315±0.170 

66±20 

0.414±0.170 



4 

1.345±0.248 

64±23 

0.425±0.248 

0.44±0.05 

J6 

1 

1.954±0.171 

71 ±20 

0.461±0.171 



2 

2.038±0.461 

59±25 

0.580±0.461 



3 

2.058±0.438 

64±27 

0.566±0.438 



4 

1.943±0.209 

53±17 

0.470±0.209 

0.00±0.10 

J7 

1 

2.641±0.152 

57±16 

0.436±0.152 



2 

2.851±0.270 

71 ±20 

0.395±0.270 



3 

2.923±0.127 

89±23 

0.414±0.127 



4 

2.990±0.176 

56±17 

0.379±0.176 

0.67±0.07 


3.2. Transverse structure 

For analyzing the transverse structure of the jets in Cygnus A, 
the four CLEAN images were convolved with a common cir¬ 
cular beam of 0.15 mas, which is approximately intermediate 
between the major and minor axes of the natural beams (Tab. 
1). Therefore, a slight super-resolution in the transverse direc¬ 
tion was applied. The analysis was performed both on the single 
epoch maps and on a stacked image (Fig. 2), created after align¬ 
ing the four maps as described in Sect. 3.1. In order to facilitate 
the study, the images were first rotated by 16 degrees clockwise, 
so that the mean jet axis was parallel to the x axis, and sliced 
pixel-by-pixel (1 px=0.03 mas) using the AIPS task SLICE in 
the transverse direction. Then, the transverse intensity distribu¬ 
tions of each slice were fitted with Gaussian profiles using the 
task SLFIT. This task was ran twice for the entire length of the 


4. Results 

4.1. Kinematic anaiysis and components iight-curves 

The proper motions p of the features described in Fig. 1 were 
calculated by fitting their radial separations z (Tab. 2, Col. 3) 
from component N with time f, as shown in Fig. 3. The motion 
of components J2 and J3 is best described by 2nd order polyno¬ 
mials, indicating an acceleration. An estimate of the kinematic 
parameters was obtained by fitting a simple parabola of the form 
z(f) = afi + c. The motion of the remaining features, instead, ap¬ 
pears more uniform and a linear fit was performed. The apparent 
speedjSapp = Vapp/c can be calculated as/3app = {pD]S)l{c{l +Z)), 
where Dl the luminosity distance, c is the speed of light and Z 
is the redshift. In the case of the accelerating features, Vapp is the 
speed calculated for the last epoch, November 2009. The values 
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2.7e-06 



Fig. 2. Top - Stacked 43 GHz image from observations in 2007-2009, rotated by 16° clockwise. It is convolved with a circular beam of 0.15 
mas. Contours represent isophotes at 0.3, 0.6, 1.2, 2.4, 4.8, 9.6, 19.2, 38.4, 76.8, 153.6 mjy/beam. The four maps were aligned to the position of 
component N, as shown in Fig. 1. Bottom - Double peaked intensity profiles at the location of the cuts indicated in the top image. 


obtained for the apparent speed of each component are reported 
in Tab. 2, Col. 6. In the approaching jet we observe an accelera¬ 
tion up toySapp = 1.15+0.04 within ~ 0.8 mas from the core. Then 
the speed drops significantly, and increases again in the outer- 
jet, but remaining largely subluminal. Three components do not 
show a significant proper motion: Cl (Jiapp = 0.00 + 0.06), JO 
(jSapp = -0.03 + 0.03) and J6 (ySapp = 0.00 ± 0.10). Given the fact 
that the entire nuclear region, represented by component N and 
by the adjacent JO and Cl, appears stationary, a slightly differ¬ 
ent choice of the reference point for the kinematic study would 
not change the results; therefore the exact core registration is in 
this case not crucial. The hypothesis that the inner-jet compo¬ 
nents are crossing an acceleration region is supported by their 
light-curves (Fig. 4). The flux densities of Jl, J2 and J3 are de¬ 
creasing with time, while those of the outer features (J4-J7) are 
approximately constant, or increasing in some cases. As the jet 
expands, a decrease of the energy emitted by the traveling fea¬ 
tures may be expected because of adiabatic losses. However, we 
argue that this effect has little impact in our case. The distance 
traveled by each component during the period of the monitoring 
is in fact small enough that they do not experience significant 
changes of the jet width. As shown in Tab. 2 (Col. 5), the sizes 
of the components are either approximately constant in time or 
they vary slightly without showing a clear trend. Instead, dif¬ 


ferential Doppler (de)-boosting may better explain the observed 
properties. It is important to keep in mind that the Doppler fac¬ 
tor (5 is a monotonically increasing function of speed /3 only for 
a tiny range of viewing angles, approximately between 0 and 
2 degrees. At larger angles, higher speeds often imply a lower 
boosting, or even a de-boosting {6 < 1). The latter means that 
the narrow boosting cone of the emission can point away from 
the observer not only in the counter-jet, but also in the approach¬ 
ing side, at a sufficiently high speed exceeding the critical speed 


P6=\ - 


2cos0 
1 + cos^ 6 


( 1 ) 


Flows featuring Lorentz factors as high as 50 are already de- 
boosted if oriented at an angle of > 11°, while a Lorentz factor 
of ~ 1.1 (y8~0.49) is sufficient at an angle of 75°. Since the view¬ 
ing angle of Cygnus A is certainly large, we find it plausible 
that the quick and strong decrease in flux density of the inner-jet 
components is rather due to a decrease of their Doppler factors, 
which strengthens the scenario suggested by the kinematic anal¬ 
ysis. In this framework, the detection of the new feature Jl in 
March 2009 gives us the opportunity to calculate the Doppler 
factor of the flow at the position occupied by Jl in November 
2009, z = 0.36 mas, and therefore to obtain an estimate for the 
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Time t [years] 


Fig. 3. Radial distance from core z versus time t for the components in 
the outer-jet (top) and those in the inner-jet (bottom). 



Time t [years] 

Fig. 4. Light-curves of the jet components. The accelerating inner-jet 
components (dashed-magenta lines) are characterized by a decreasing 
flux density with time, while the outer ones (continuous lines-other col¬ 
ors), moving at constant speed, feature a more constant flux. 


jet viewing angle 0. Assuming that the flux density of J1 is close 
to the intrinsic one in March 2009 (due to the low intrinsic speed 
we expect at that position) the decrease of flux density constrains 
6 according to: 



where a is the spectral index, 5obs is the flux density of J1 in 
November 2009 and 5int is that in March 2009. With a classical 
synchrotron spectral index a - -0.7 {S oc v"), we obtain 6 - 


0.92. A second constraint on the apparent speed jS^pp = 0.75 of J1 
in November 2009 allows us to solve a system of 2 equations for 
deriving the viewing angle 6 and the intrinsic speed j3 at z=0.36 
mas: 


j8sin6> 

1-/JCOS0’ l-yScos0 


(3) 


obtaining [5{z - 0.36) = 0.65 and 6 - 74.5°. The value of the 
viewing angle obtained can be considered an upper limit, be¬ 
cause if the change of Doppler factor is not the only cause for 
the decrease of flux, a smaller viewing angle is required. 


4.2. Ridge-line 

The pixel-based analysis of the transverse intensity profiles 
shows that they are double peaked for most of the length of the 
source (Fig. 5). The limb brightening is most pronounced in the 
inner-jet, while a single ridge is visible further out, at a separa¬ 
tion from the core between 2 and 4 mas. In the counter-jet we 
also find evidence for a limb brightened structure: the two ridge 
lines are closer compared to the main jet, but still well separated 
after accounting for positional errors, assumed equal to one fifth 
of each FWHM. An error bar which also takes account of the 
SNR of the profile (as described in Fomalont 1989) was calcu¬ 
lated for comparison. The latter yields unreasonably small errors 
in the brightest regions, while in the faint ones they are anyhow 
smaller compared to the first method of calculation. Therefore 
the most conservative error bars were chosen. In the central ~ 
0.5 mas, the poor resolution does not allow a double profile to 
be clearly traced, and results from a single Gaussian fit are plot¬ 
ted. The same is done in the other regions of the jet where a 
well defined double ridge is missing (error bars of the two peaks 
overlapping or bad fit). The disappearance of the double pro¬ 
file in the main jet (at z~ 2 mas) is accompanied by an apparent 
bending of the jet axis towards the south and by a narrowing of 
the jet width, with the two ridges smoothly approaching each 
other. Interestingly, this position coincides with the location of 
the stationary component J6. While additional interpretations of 
this feature are proposed in Sect. 5.3, a more detailed analysis is 
necessary to establish if the observed structure can be ascribed to 
helical Kelvin-Helmoltz or to current-driven instabilities devel¬ 
oping in the jet. However, following the previous discussion on 
relativistic effects in flows seen at large viewing angle (in sum¬ 
mary: “the faster, the dimmer”), the most natural explanation for 
the observed limb brightening is the presence of a transverse gra¬ 
dient of the bulk Lorentz factor of a spine-sheath kind, i.e. with 
F decreasing from the jet axis towards the edges. 


5. Discussion 

5.1. Two-layers acceleration 

The kinematic and light-curves analysis of the VLBI maps at 7 
mm has revealed the existence of an acceleration region in the 
inner-jet extending up to a de-projected distance of ~0.9 pc, as¬ 
suming the viewing angle 0 = 74.5° derived in Sect. 4.1. In the 
outer-jet, lower speeds are determined, in agreement with find¬ 
ings from VLBI at lower frequencies (Bach et al. 2003). The ob¬ 
served fading of the fast flow (Fig. 4) suggests that these lower 
speeds do not represent an intrinsic deceleration, but are sim¬ 
ply the signature that the contribution to the emission from the 
slower external layers has become dominant in that region. If 
we identify the two regions as representing the “fast” and the 
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Fig. 5. Ridge line structure of Cygnus A at 7 mm. The magenta color is used in those regions were a well defined double ridge line is present, while 
green points describe the remaining single-peaked regions. The error bars on the peak positions are set to one fifth of the de-convolved FWHM. 



“slow” layers of the flow respectively, we can study separately 
their kinematic properties. In Fig. 6 (top) we plot the calcu¬ 
lated intrinsic speeds /3 versus the de-projected distance from the 
core expressed in units of Schwarzschild radii Rs. We then note 
that the slow flow is also accelerating, but with a milder gra¬ 
dient compared to the fast one. In this case though, the range 
of speeds is such that the Doppler factor remains approximately 
constant, so we do not observe strong variations of the flux den¬ 
sity. Both sections are accelerating over a large distance, at least 
of the order of ~ 10^ - 10'*/?s. This is a characteristic signature of 
relativistic magnetically driven outflows (Begelman & Li 1994; 
Contopoulos 1995; Vlahakis & Konigl 2004), in which the con¬ 
version from Poynting to kinetic flux is found to proceed quite 
slowly. The steeper gradient of speed which we measure for the 
fast flow (central spine) compared to the slow one (sheath) in¬ 
dicates that the conversion to kinetic energy proceeds faster for 
sections of the flow which are closer to the jet axis. By fitting a 
power law (Fig. 6, top) we And oc 20-5i±o.04 former, and 

P oc 20-43 ±o.oi jjjg jatter. This translates into a quite different 
behavior of the Lorentz factor (Fig. 6, bottom), whose depen¬ 
dence on the core separation is F oc 20-50±o.i9 p ^o.ii±o.o4 
for the two flows respectively. 


5.2. Mean opening angle 

From our study, other important parameters of the flow, like the 
shape and opening angle, can be estimated. From the analysis 
of the transverse structure (Sect. 3.2), we can infer the depen¬ 
dence of the jet width on the core separation. It is preferable in 
this case to take into account not the stacked map but the sin¬ 
gle images, as this prevents the temporal evolution of the ridge 
line from affecting the results. Here we present our findings for 
the last epoch, November 2009, which is the one showing most 
prominently the double ridge line structure. The jet width was 
calculated for each slice as the separation between the edges of 
the two (de-convolved) FWHM. The results are shown in Fig. 7, 
where the de-convolved widths obtained from a single-Gaussian 
slice fit are plotted for comparison, and show good agreement. 
This indicates that the expansion of the two ridge lines coin¬ 
cides well with the overall expansion of the jet. We also plot the 
widths obtained from the modelflt analysis of all the four ob¬ 
serving epochs. With some exceptions (discussed below) there 
is again a substantial overlap. The jet width is clearly increas¬ 
ing both in jet and counter-jet up to approximately 2 mas from 
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Fig. 6. Top - Intrinsic speed as a function of de-projected core separa¬ 
tion. A fast and a slow component with different acceleration gradients 
can be identified. Bottom - Corresponding Lorentz factor versus de- 
projected core-separation. 

component N (at z = 0 in the plot). By performing a linear fit 
in this region (|z| < 2) we can calculate the mean apparent full 
opening angle 0appj increased relative to the intrinsic one 0 by a 
factor 1/ sin 6 due to projection. In jet and counter-jet we obtain 
respectively; 

= 10.2° + 0.3°, = 4.8° + 0.4° (4) 

which for 6 = 74.5° correspond to intrinsic opening angles; 

<^j = 9.8° + 0.3°, = 4.7° + 0.4° (5) 

As found also in Krichbaum et al. (1998) at 22 GHz, the jet ap¬ 
pears about two times broader than the counter-jet. Following the 
classical assumption of intrinsic symmetry, the apparent open¬ 
ing angle should in principle be the same in jet and counter-jet, 
unless they are strongly misaligned. Although a small misalign¬ 
ment is present both on pc and kpc scales, this is not enough 


Article number, page 6 of 9 











B. Boccardi* et al.: The stratified two-sided jet of Cygnus A 


to justify the observed difference. A possible explanation could 
be the different SNRs of the approaching and receding sides, 
with the counter-jet being considerably dimmer than the jet at 
the same core separation and thus detectable over a narrower 
section. Also, we note that there is, in this case, a certain dis¬ 
crepancy between the modelfit size of Cl and the correspond¬ 
ing width of the intensity profile, with the first being larger and 
closer to the width of the jet side. This is probably an effect 
from the smaller beam used in the slicing and it suggests that 
the counter-jet emission is less compact than the jet emission. 
However the hypothesis of intrinsic asymmetry, indeed possible 
(e.g Fendt & SheikJinezami 2013), cannot be excluded. 

Based on these results, the jets in Cygnus A can anyhow be 
defined relatively broad on the examined scale of ~2 mas. It is 
useful to note that the fast jet expansion cannot be clearly per¬ 
ceived by eye in the map in Fig. 2. This is first of all due to 
the beam convolution effect, causing especially the base of the 
jet to appear much wider than it actually is. The second rea¬ 
son is that, even after de-convolving the jet width, an apparent 
full opening angle of “only” ~10° is difficult to appreciate visu¬ 
ally on the scale of the image. Despite the misleading appear¬ 
ance, however, the jets in Cygnus A are considerably broader 
than in blazars, in which intrinsic opening angles of ~(1 - 2)° 
are typically measured (Pushkarev et al. 2009). This result, also 
found by examining other radio galaxies on parsec-scale, repre¬ 
sents another piece of evidence for the spine-sheath scenario, in 
which sources seen at larger viewing angles become more and 
more sheath-dominated, and therefore broader, while the emis¬ 
sion from blazars mostly comes from the strongly boosted, thin 
spine. 


De-projected z [Rg] 
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Fig. 7. Jet aperture versus core separation in November 2009 from the 
double ridge line analysis (magenta), from a single Gaussian fit (green). 
Overplotted in blue are the sizes of the modelfit components from all 
epochs. Component J2 occupies only a portion of the jet in the trans¬ 
verse direction. Note that, in order to more clearly show the expansion, 
the X and y axis have different scales. 


5.3. Narrowing of the jet at z~2 mas 

The increase of the jet width is not always smooth, but is charac¬ 
terized by some oscillations. Although an oscillating expansion 
can occur in jets (Daly & Marscher 1988), we cannot exclude 
that this is an image artifact, especially in the faint regions of 
the counter-jet. The narrowing at ~ 2 mas in the main jet, and 
perhaps in the counter-jet too is, however, very prominent and 
coincides with the location of the stationary feature J6. Station¬ 
ary features in jets can be related to instability patterns, as men¬ 


tioned in Sect. 4.2, or may result from re-collimation mecha¬ 
nisms. In the purely hydrodynamic case, recollimations naturally 
occur when the expanding flow becomes under-pressured com¬ 
pared to the medium (Daly & Marscher 1988). In MHD flows, 
collimation features can also appear when the jet ceases to be 
strongly magnetized and causally connected to the central en¬ 
gine (Meier 2012). The transition to the new regime is possibly 
accompanied by a saturation of the acceleration, which is nev¬ 
ertheless not observed in this region. Kinematic VLBI studies 
at lower frequencies (Bach et al. 2003, and references therein) 
have shown that the acceleration in the main jet extends at least 
to ~5 mas, where another prominent narrowing appears (Krich- 
baum et al. 1998). The nature of this collimation feature is ulti¬ 
mately not clear. It is interesting though that the narrowing is not 
present if we look at the modelfit widths (Fig. 7). In other words, 
the larger beam used in the modelfit analysis does not “see” the 
feature, indicating that whatever process is taking place, it must 
involve the flow on smaller scales. 

5 . 4 . Collimation regime 

Theoretical models and simulations describing the conversion 
from Poynting flux to kinetic flux in magnetically-driven jets 
have shown that the characteristics of the external medium 
represent a fundamental parameter (Komissarov et al. 2007; 
Tchekhovskoy et al. 2008; Lyubarsky 2009). Specifically, the 
gradient of the external pressure along the jet largely deter¬ 
mines the shape of the jet and the acceleration rate, which are 
intimately connected. In an analytical study of the asymptotic 
structure of magnetically-driven flows, Lyubarsky (2009) pre¬ 
dicts the existence of two main regimes. In the first, the equilib¬ 
rium regime, the pressure of the external medium p decreases 
slowly, as p oc with k < 2. In these conditions the magnetic 
field in the flow succeeds in keeping a strong poloidal compo¬ 
nent. The jet radius r as a function of distance from the core 
is then given by a power law with index smaller than 0.5, the 
Lorentz factor increases as the radius (F oc r) and the asymptotic 
£ shape of the jet is cylindrical. If the pressure gradient is steeper, 
1 with index k > 2, the poloidal magnetic field becomes negli- 
J gible, while the toroidal component cannot provide strong con¬ 
finement to the flow, which then gets broader, accelerates more 
slowly (F oc z/r), and eventually reaches a conical shape. This 
is the non-equilibrium regime. The transition state between the 
two regimes occurs for k = 2. In this case, the asymptotic be¬ 
havior of the flow is determined by the intensity of the magnetic 
field Bq and by the external pressure po at the base of the flow. 
Specifically, if the ratio (hnpojlB^ is larger than 0.25, the flow 
has characteristics similar to the first regime, and has a perfectly 
parabolic shape, i.e. with a radius r increasing as A smaller 
ratio instead leads the flow to a state similar to the second regime. 

In the following, we show that Cygnus A resembles the prop¬ 
erties expected in the transition state {p oc z~^) and that the 
asymptotic regime is that of equilibrium. Let us consider the 
shape of the jet. Because of the higher SNR, we focus our anal¬ 
ysis on the approaching side. In order to reduce the impact of 
small scale oscillations and to study the overall expansion of the 
jet, we consider the modelfit sizes in the four epochs up to the 
narrowing at ~ 2 mas. We exclude component J2 because, as 
shown in Fig. 7, this is clearly a feature occupying, transver- 
sally, only a portion of the jet. Results are shown in Fig. 8. The 
expansion can be described by a single power law of the form 
r oc 20-55 ±o.07 jjp J.Q i-C- the jet has a parabolic shape. 

Two elements suggest that the flow collimates in the equilibrium 
regime. Firstly, the Lorentz factor gradient of the fast component 
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calculated in Sect. 5.1 (F cx 20-50±o.i9^ compatible with the pre¬ 
diction that r grows as r, at least in the central body of the jet. 
Secondly, VLBI studies at 5 GHz (Carilli et al. 1991) showed 
that the jet has a cylindrical shape between 4 and 20 mas, which 
is the asymptotic state in the equilibrium regime. Interestingly, 
4-5 mas (~2.5 x 10'*/?s) is also the maximum scale on which 
acceleration has been observed in this source (Bach et al. 2003), 
and where a second, prominent recollimation occurs (Krichbaum 
et al. 1998). Therefore this location can be identified as a good 
candidate for a region where the jet of Cygnus A ceases to be 
strongly magnetized and becomes causally disconnected from 
the central engine. 

De-projected z [Rg] 

1000 10000 



cc 


T 1000 


the most narrow region of the jet. Instead, both jet and counter¬ 
jet appear to emanate exactly from the gap of emission, which 
features the absolute minimum width. It is then likely that the 
gap marks the true location of the central engine, which would 
then be located ~800 Rs upstream of the jet base, represented 
by component N. Given the amount of evidence for the jet being 
magnetically driven, it is interesting to wonder how the jet would 
appear when it is still completely magnetically dominated, i.e. 
with a magnetization parameter much larger than 1. Due to the 
very low speed the flow would have, it would also be faint and 
may appear as a gap. Other explanations are of course possible. 
The gap could alternatively be the effect of strong absorption by 
a molecular torus, which is certainly present and obscures part 
of the counter-jet (Struve & Conway 2010). Another possibility 
is that this is simply a self-absorbed region. A detailed spectral 
analysis using higher frequency data is required and in progress. 


_ 6. Conclusions 


De-projected z [mas] 

Fig. 8. Power-law fit of the jet width (FWHM=2*radius r) versus de- 
projected distance from the core, as inferred from MODELFIT. 


5.5. Comparison with M87 

In the previous section we have shown that the prototypical FRII 
radio galaxy Cygnus A is collimated and accelerated on a scale 
of ~2.5 X lO^Rs^ with the jet being parabolic in this region and 
cylindrical further downstream. It is interesting to compare these 
properties with those inferred for M87, the most studied FRI ob¬ 
ject. Asada & Nakamura (2012) have shown that the jet of M87 
is also parabolic up to a scale of ~ lO^Rs, with r oc thus with 
a power law index very close to that we measure for Cygnus 
A. However, after the HST-1 feature, the jet becomes conical 
rather than cylindrical, just as expected in the non-equilibrium 
regime. An important difference between the equilibrium and 
non-equilibrium regimes resides in the efficiency of conversion 
of Poynting flux. While in the former case the Poynting flux 
will be mostly converted into kinetic flux, until equipartition 
is reached, a jet in non-equilibrium will remain magnetized up 
to large distances without reaching very high Lorentz factors. 
The high power, thin jets of the FRII sources and the lower 
power, broader jets of the FRIs appear to fit well into the two 
categories respectively, indicating that the FRI/FRII dichotomy 
might result from different environmental conditions. The dis¬ 
covery of sources with hybrid morphologies (Gopal-Krishna & 
Wiita 2002), the so called HYMORS, supports the same sce¬ 
nario. 


5.6. Nature of the gap of emission 

Finally, we comment on the nature of the sharp gap of emis¬ 
sion situated at ~ 0.2 mas east of component N and of ~ 0.1 
mas in size. Figure 7 shows that, although the zero in the x axis 
is set to coincide with the location of component N, this is not 


We have performed a high resolution VLBI study at 7 mm of the 
kinematic properties and transverse structure of the two-sided jet 
in Cygnus A. Results can be summarized as follows. 

- We detect an acceleration of the flow up to supeiiuminal 
speeds (jSU^“ = 1.15 ± 0.04) in the first parsec of the jet. 
The outer jet exhibits lower speeds, in agreement with stud¬ 
ies at lower frequencies. We suggest that a fast and a slow 
component exist in the flow, and interpret the slower speeds 
in the outer jet as an effect from the de-boosting of the fast 
component, naturally occurring at large viewing angle. 

- The acceleration region of both the fast and slow component 
is extended (~10^ - lO^Rs); as characteristically predicted 
for magnetically-driven outflows. 

- Both jet and counter-jet are limb brightened. In the presence 
of a transverse gradient of the bulk Lorentz factor, the limb 
brightening is evidence for a spine-sheath structure of the jet, 
as any fast flow appears faint at large viewing angle. 

- The mean opening angles in the two sides are large (the full 
angle in the jet is ~10°) compared to typical values inferred 
for blazars (1° -2°). This can be explained by the possibility, 
in radio galaxies, to observe not only the thin inner spine but 
also the outer sheath. 


- The shape of the flow is parabolic (r oc z* 


, 0 . 55 ± 0.07 


). The ac¬ 


celeration gradients and the collimation regime can be well 
reproduced by assuming an external pressure p which scales 
with distance from the central engine z as p oc and the 
ratio (67r/7o)/Bo to be larger than 0.25 at the base of the flow, 
as described in Lyubarsky (2009). 

- The jet stops being strongly magnetized possibly on a scale 
of ~2.5 X 10"^Rs^ where it enters the asymptotic cylindrical 
shape (Carilli et al. 1991), the bulk of the acceleration sat¬ 
urates (Bach et al. 2003) and a recollimation occurs (Krich¬ 
baum et al. 1998). 
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